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Abstract. Using the methods originally developed for Random Matrix Theory we derive an exact mathe- 
matical formula for number variance A«(L) (introduced in [4]) describing a rigidity of particle ensembles 
with power-law repulsion. The resulting relation is consequently compared with the relevant statistics of 
the single-vehicle data measured on the Dutch freeway A9. The detected value of an inverse temperature 
P, which can be identified as a coefficient of a mental strain of car drivers, is then discussed in detail with 
the respect to the traffic density q and flow J. 



PACS. 05.40.-a Fluctuation phenomena, random pi 
portation - 05. 45. -a Nonlinear dynamics and chaos 

1 Terminus a quo 

The main goal of this paper is to show that the statistical 
fluctuations of single-vehicle data (in vehicular flows) can 
be very well predicted by the methods known from Ran- 
dom Matrix Theory where used for describing the statis- 



:esses, noise, and Brownian motion - 89.40.-a trans- 
tics of energy levels in quantum chaotic systems. Above 
that, we intend to demonstrate that the changes of statis- 
tical variances in vehicular samples depends not only on 
macroscopical traffic quantities and three traffic phases (as 
published in p] or [2]), but also on psychological character- 
istics of driver's decision-making process. We will evince 
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that transitions among the traffic phases (free flow, syn- and 
chronized flow, and wide moving jam) cause perceptible 
changes in the mental strain of car drivers. 



r p(r) d?' = 1. (5) 



The latter represents a scaling to the mean clearance equal 
As reported in Ref. [I], 0, [3], g] and [5] three traffic to on6j which ig introduced for convenience. The above- 
phases show substantially different microscopical proper- men ti ncd distribution © is in a good agreement with 
ties. It was demonstrated in [6], [7], [8], and £U\ that the thc clearance distribution of real-road data (0, [7], g], 
microscopical traffic structure can be estimated with the and jgj) wnC reas thc inverse temperature is related to 
help of a repulsive potential (applied locally only) describ- thc traffic dcnsity Q . Wc notc that thc inverse temperature 
ing thc mutual interaction between successive cars in the p can be unders tood as a quantitative description of thc 
chain of vehicles. Especially, thc probability density p(r) mental strain of drivcrs in a givcn traffic situation. More 
for thc distance r of two subsequent cars {clearance dis- spec ially, the parameter reflects thc psychological pres- 
tribution) can be derived by means of an one-dimensional sure levcl under which the driver is if mov ing in thc traffic 
gas having an inverse temperature and interacting by strcam . whereas thc free flows induces practically no mcn- 
the repulsive potential V(r) = r" 1 (as discussed in Ref. tal strain of driverS; with tnc increasing traffic density thc 
0, [6], and [8]). Concretely, denoting &{x) the Heaviside's mental strain esca lates rapidly. This will be confirmed in 

function the following text. 

I I, x > 
0(x) = \ 

[o, i^O, 

and JC\{x) the modified Bessel's function of the second 
kind of order A (Mac-Donald's function), the clearance dis- 
tribution of the above-mentioned thermodynamical traffic 
gas reads as 



(r) =AO(r)e-h- Br , (1) 



2 Statistical variances in traffic data 

Another powerful way to inspect thc interactions between 
cars within the highway data is to investigate the traffic 
flow fluctuations. One possibility is to use the so-called 



where time-gap variance Z\ T considered in paper [3] and defined 

3 _ e -VP 

B — H , (2) as follows. Let {ti : i = 1 . . . Q} be the data set of time 

, , 10 , i , , „ intervals between subsequent cars passing a fixed point on 

<AT 1 = 2W^/Ci(2 v /B/3). (3) vox 

thc highway. Using it one can calculate thc moving average 

We remark that p(r) fulfils two normalization conditions 



T (N) 



p(r)dr = l (4) 1 k ~ N 



k+N-l 

- Y, ^ (k = l...Q-N+l) 
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of the time intervals produced by the N + 1 successive 
vehicles (i.e. N gaps) as well as the total average 



T 



Q 



i=l 



The time-gap variance A T is defined by the variance of 
the samplc-avcragcd time intervals as a function of the 
sampling size N, 

Q-N+l 



Aj 



Q-N 



I V (t {n) 

n + i v k 

fe=l 



T 



where k runs over all possible samples of N + 1 successive 
cars. For time intervals ti being statistically independent 
the law of large numbers gives A T (N) ex. 1/N. 

A statistical analysis of the data set recorded on the 
Dutch freeway A9 and published in Ref. [9] leads, how- 
ever, to different results - see the Figure 1. For the free 
traffic flow (g < 15 veh/km/lane) one observes indeed the 
expected behavior Aj(N) oc 1/N. More interesting be- 
havior, nevertheless, is detected for higher densities (g > 
35 veh/km/lanc). Here Nishinari, Treiber, and Hclbing (in 
Ref. [5]) have empirically found a power law dependence 

A T (N) cx N e 

with an exponent e f=a —2/3, which can be explained as a 
manifestation of correlations between the queued vehicles 
in a congested traffic flow. 

There is, however, one substantial drawback of this 
description. The time-gap variance was introduced ad hoc 
and hardly anything is known about its exact mathemati- 
cal properties in the case of interacting vehicles. It is there- 
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Fig. 1. The time-gap variance A T (N) as a function of the sam- 
pling size N (in log- log scale). Plus signs and stars represent 
the variance of average time-gaps for free and congested flows, 
respectively. 

fore appropriate to look for an alternative that is math- 
ematically well understood. A natural candidate is the 
number variance A^(L) that was originally introduced for 
describing a spectral rigidity of energy levels in quantum 
chaotic systems, i.e. for describing a structure of eigenval- 
ues in the Random Matrix Theory. An(L) reproduces also 
the variances in the particle positions of a certain class of 
one-dimensional interacting gases (for example Dyson gas 
in Ref. [12) and it is defined as follows. 

Consider a set {r^ : i = 1 . . . Q} of distances (i.e. clear- 
ances in the traffic terminology) between each pair of sub- 
sequent cars moving in the same lane. We suppose that 
the mean distance taken over the complete set is re-scaled 
to one, i.e. 

Q 

]>> = Q. 

i=l 

Dividing the interval [0, Q] into subintervals [{k— 1)L, kL] 
of a length L and denoting by nk(L) the number of cars 
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in the kth subinterval, the average value n(L) taken over 

all possible subintervals is 

1 IQ/H 

= IqTZJ g ML) = u 

where the integer part \Q / L\ stands for the number of all 
subintervals \(k — l)L,kL] included in the interval [0,(5]. 
Number variance As(L) is then defined as 

1 LQ/iJ 
^) = IQ7ZJ g { n k{L) - L) 2 

and represents the statistical variance of the number of 
vehicles moving at the same time inside a fixed part of 
the road of a length L. The mathematical properties of 
the number variance are well understood and therefore 
A N (L) serves as a better alternative for description of traf- 
fic fluctuations (a rigidity of the vehicular chain) than the 
time-gap variance A T (N) itself (see also [1]). 

3 Exact formula for number variance 

Considering the probability density ([TJ) with the only free 
parameter (3 (the inverse temperature) we aim to derive an 
exact formula for the rigidity A$(L) of thermodynamical 
traffic gas (see [7]). For these purposes we use the methods 
presented in [12] and [11] in detail. 
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Fig. 2. The slope x(P) °f the number variance Aa(L) = 
as a function of the inverse temperature j3. The squares rep- 
resent the results of numerical calculations whereas the con- 
tinuous curve displays the result ([8} obtained by the exact 
mathematical computations. The behavior close to the origin 
is magnified in the inset. 

calculated via recurrent formula 

Pn(r) = p n -i(r) * po(r), 

where symbol * represents a convolution of two indepen- 
dent probabilities, i.e. 

Pn(r)= / p n -i(s)po(r - s)ds. 

Using the approximation of the function 



9n(s) = e 



Let p n (r) represent the n-th nearest- neighbor prob- in the saddle point one can obtain 
ability density (n-th probability density for short), i.e. 
the probability density for the spacing r of n + 2 neigh- 
boring particles. Owing to the fact that probability den- where 
sity for spacing between two succeeding particles (cars) is 
p(r) = po(r), the n-th probability density p n {f) can be 



-Br 



n+1 



AC 1 = 2 



(n + 1) K n+1 (2{n + l)y/B0) 
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fixes up the proper normalization J R p n (f) dr = 1. In ad- 
dition to that the mean spacing equals 



0.2 



rp n (r) dr = n + 1. 

According the book [12] the variance A N (L) could be eval- 
uated by the formula 



A„(L) = L — 2 (L- r)(l - R(r)) dr, 



(6) 



where 



fl(r) = ^>„(r) 

n=0 

is the two-point cluster function. The convenient way to 
calculate the asymptotic behavior of the variance A^(L) 
for large L is the application of the Laplace transformation 
to the two-point cluster function y(t) — L R(r) e~ rt dr. It 
leads to a partial result 

B \^ K„ +1 (2(«+l)#+f) 



y(t) = £ 



n=0 



B + t 



The small— a; asymptotic behavior of the Mac-Donald's 
function 

2 n-X ; 



l r(n) x 

■ e 



(x < 1) 



1 



B e 2 SW 

Applying the Maclaurin's expansion (Taylor's expansion 
about the point zero) of the function h(t) = t ■ y(t) to 
order t 2 we obtain 

1 



y(t) w j + ao + c^t + Oit 2 ), 



where 



a 



2Bf3 + 3y/Bj3 
4(l + v r B^) 2 ' 
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Fig. 3. The shift 7(/3) of the number variance zAn(£) = X-^ + 7 
as a function of the inverse temperature /3. The squares rep- 
resent the results of numerical calculations whereas the con- 
tinuous curve visualizes the result (J9]) obtained by the exact 
mathematical computations. The behavior close to the origin 
is magnified in the inset. The dash-dotted line displays the 
asymptotic tendency in 7 = j(/3), i.e. Iirri0_, oo j(/3) = 1/6. 



01 



6y/B(3 + B(3(2l + 45/3 + 16y/E0) 



48S(l + 2VS/3) 3 
Then we get from equation © 



A 1I (L)= X L + J + 0(L- 1 ), 



(where .T(x) represents the gamma-function) provides the where 
approximation 

' B + t e V(^^ 



X = X(P) 



fB[3 



2B(l + VBP) 



(7) 



(8) 



and 



7 = 7 (/3) = 



6/8/3 + 5/3(21 + 45/3 + 16/5/3) 



(9) 



24(1 + 

This crowns the effort to derive an exact form for the 
number variance. The linear dependence represents 
a large— L approximation and its correctness is demon- 
strated in the Fig. 2 and Fig. 3 where compared to the 
results of numerical computations. 
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4 Number variance of traffic data 



As already discussed the behavior of the number variance 
is sensitive to the temperature (3 - or in the terminology 
of the Random Matrix Theory - to the universality class 
of the random matrix ensemble. To use the known math- 
ematical results one has not to mix together states with 
different densities - a procedure known as data unfolding 
in the Random Matrix Theory. For the transportation this 
means than one cannot mix together traffic states with dif- 
ferent traffic densities (as published in and hence with 
a different vigilance of the drivers. So we will perform a 
separate analysis of the data-samples lying within short 
density intervals to prevent so the undesirable mixing of 
the different states. 
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Fig. 4. The number variance A S (L) as a function of the length 
L. Plus signs, diamonds, stars, circles, crosses, and squares 
represent the number variance of real-road data in the selected 
density regions (see the legend for details). The curves show 
the linear approximations of the individual data. Their slopes 
were carefully analyzed and consecutively visualized in the Fig. 
5 (top part). 



Anyway, we divide the region of the measured densities 
g € [0, 85 veh/km/lane] into eighty five equidistant subin- 
tervals and analyze the data from each one of them sep- 
arately. The number variance A$(L) evaluated with the 
data in a fixed density interval has a characteristic lin- 
ear tail (see Fig. 4) that is well known from the Random 
Matrix Theory. Similarly, such a behavior was found in 
models of one-dimensional thermodynamical gases with 
the nearest-neighbor repulsion among the particles (see 
Ref. We remind that for the case where the interac- 

tion is not restricted to the nearest neighbors but includes 
all particles the number variance has typically a logarith- 
mical tail - see [T2]. So the linear tail of A$(L) supports 
the view that in the traffic stream the interactions are re- 
stricted to the few nearest cars only. The slope of the linear 



tail of A N (L) decreases with the traffic density (see the top 
subplot in the Fig. 5). It is a consequence of the increas- 
ing alertness of the drivers and hence of the increasing 
coupling between the neighboring cars in the dense traffic 
flows. 

The exactly derived properties of the function A N (L) 
agree with the behavior of the number variance extracted 
from the traffic data (see Fig. 4). A comparison between 
traffic data number variance and the formula |(7J) allows 
us to determine the empirical dependence of inverse tem- 
perature P on traffic density g. The inverse temperature 
reflects the microscopic status in which the individual ve- 
hicular interactions influence the traffic. Conversely, in the 
macroscopic approach, traffic is treated as a continuum 
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Fig. 5. The slope \ and the inverse temperature j3 as func- 
tions of the traffic density g. The squares on the upper sub- 
plot display the slope of the number variance Au(L) (see 
Fig. 4), separately analyzed for various traffic densities. The 
lower subplot visualizes the fitted values of the inverse tem- 
perature f3, for which the exact form of number variance 
A S (L) = x(fl) L + 7(/3) corresponds to the number variance 
obtained from the traffic data. The continuous curves repre- 
sent polynomial approximations of the relevant data. 

and modelled by aggregated, fluid-like quantities, such as 
density and flow (see pQ and [5] ) ■ Its most prominent result 
is the dependence of the traffic flow on the traffic density 
- the fundamental diagram. 

It is clear that the macroscopic traffic characteristics 
are determined by its microscopic status. Consequently 
there should be a relation between the behavior of the 
fundamental diagram and that of the inverse temperature 
(3. In the Figure 6 we display the behavior of the inverse 
temperature (3 simultaneously with the fundamental dia- 
gram. The both curves show a virtually linear increase in 
the region of a free traffic (up to g w 10 veh/km/lane). 



The inverse temperature (3 then displays a plateau for the 
densities up to 18 veh/km/lane while the flow continues 
to increase. A detailed inspection uncovers, however, that 
the increase of the traffic flow ceases to be linear and be- 
comes concave at that region. So the flow is reduced with 
respect to the outcome expected for a linear behavior - 
a manifestation of the onset of the phenomenon that fi- 
nally leads to a congested traffic. For larger densities the 
temperature (3 increases up to g 32 veh/km/lane. The 
center of this interval is localized at g w 25 - a criti- 
cal point of the fundamental diagram at which the flow 
starts to decrease. This behavior of the inverse tempera- 
ture is understandable and imposed by the fact that the 
drivers, moving quite fast in a relatively dense traffic flow, 
have to synchronize their driving with the preceding car 
(a strong interaction) and are therefore under a consid- 
erable psychological pressure. After the transition from 
the free to a congested traffic regime (between 40 and 
50 veh/km/lane), the synchronization continues to decline 
because of the decrease in the mean velocity leading to de- 
creasing (3. Finally - for densities related to the congested 
traffic - the inverse temperature increases while the flow 
remains constant. The comparison between the traffic flow 
and the inverse temperature is even more illustrative when 
the changes of the flow are taken into account. Therefore 
we evaluate the derivative of the flow 

. dJ 



J' 



dg 



The result of the evaluation can be seen from the Figure 
5 where one can trace the significant similarities between 
the shape of J' = J'(g) and inverse temperature f3. 
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ing known from Random Matrix Theory) revealed that 
the rigidity of traffic samples (evaluated by means of the 
number variance shows a linear dependence 
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Fig. 6. Traffic flow J(g) and inverse temperature /3(g) as func- 
tions of the traffic density g. Plus signs display a traffic flow 
in thousands of vehicles per hour and squares correspond to 
inverse temperature of the traffic gas. The results of a polyno- 
mial curve-fitting are visualized by the continuous curves. The 
dashed lines represent a linear approximations of the relevant 
data near the origin. 

5 Summary and conclusion 

With the help of methods of Random Matrix Theory we 
have derived the exact formula for the rigidity (quantified 
by the number variance A$(L)) of the one-dimensional 
thermal particle-gas with the repulsion-potential which 
depends on the reciprocal value of clearances among the 
nearest-neighboring particles. The correctness of the result 
obtained was confirmed by the comparison to the numer- 
ical solution. 

Analogously to the investigations presented in [4] we 
have analyzed the number variance An of single-vehicle 
data measured on the Dutch freeway A9. The robust sta- 
tistical analysis (based on the principle of data unfold- 



(10) 



(in each of eighty five equidistant density-subregions) whose 
slope x depends on traffic density significantly. This knowl- 
edge substantially specifies the general results published 
in Ref. 0. 

Above that, the comparison of the exact result (JTUJ) 
with the number variance of traffic data has provided the 
remarkable possibility how to detect the mental strain 
level of car drivers moving in the traffic stream, ft was 
demonstrated that the inverse temperature of the traffic 
sample, indicating the degree of stress of the drivers, shows 
an increase at both the low and high densities. In the in- 
termediate region, where the free flow regime converts to 
the congested traffic, it displays more complex behavior. 

To conclude, we have shown that microscopical struc- 
ture of traffic ensembles is rapidly changing with the traffic 
density and is very well described by the one-parametric 
class of distributions JT]) where the only parameter is the 
mental strain coefficient /3 depending on traffic density g. 
In addition, a recent research (Ref. [X5] ) shows that a suit- 
able modification (using the knowledge presented in this 
paper) of Nagel-Schreckenberg cellular model would lead 
to an interesting progress in the traffic modelling. 
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Fig. 7. The inverse temperature /3(g) and the first derivative 
J' = J'(q) as functions of the traffic density q. Squares corre- 
spond to the inverse temperature of traffic sample while plus 
signs display the average value of J' = J'(q) ( m kilometers 
per hour). The continuous curves represent relevant polyno- 
mial approximations. 
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